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Background: Close Formations 




14% reduction in heart rate when in formation (Welmerskirtch, et. al, 2001) 
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Over 20% drag reduction and 18% fuel flow reduction for trailing aircraft (NASA, 2002) 
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Extended Formations 


• ~15-40 span separation 


• A safer approach to formation flight 


• Can be implemented with today’s 
aircrafts without modification 



• Up to 10% fuel flow reduction for 
transport aircraft* 


*Pahle, Joe, “A preliminary Investigation of Formation Flight for Drag Reduction on the C-17 aircraft” 
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• Provides a safer approach to formation flying, yet retains much of the 
fuel savings seen in close formation flight 
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Extended Formations 



• Provides a safer approach to formation flying, yet retains much of the 
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• Provides a safer approach to formation flying, yet retains much of the 
fuel savings seen in close formation flight 
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Modeling 



• In close formation, the trailing aircraft can influence the lead 

• At 15-40 spans separation, this effect becomes negligible 

• Extended formations decouple lead and follower aircraft 



I I 
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Modeling 


• In close formation, the trailing aircraft can influence the lead 

• At 15-40 spans separation, this effect becomes negligible 

• Extended formations decouple lead and follower aircraft 



Lead Aircraft Trail Aircraft 
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Modeling Approach 
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Modeling Approach 
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Modeling Approach 
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Modeling Approach 
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Modeling Approach 





Trim 
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Out-of-Formation 
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In-Formation 




Incoming 

vortex 
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In-Formation 



Incoming 

vortex 
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In-Formation 




Roll-trim is a critical component of 
formation flying 


Incoming 

vortex 
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Roll-Trim Strategy 


( ; >| 

Trail Aircraft Euler 

Computation 


• Algorithm for trim: 


(1) Compute trail A/C flow solution 


(2) Check Lift & Roll tolerance 


(3) Deflect ailerons if necessary (roll) 

(4) Adjust AOA if necessary (lift) 

/ \ 

Rotate geometry to new 
AOA, re-mesh 

\ ) 


/ \ 

Rotate Ailerons w.r.t. 
geometry 

\ ) 


Check Loads 
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Flow Analysis: Objective Functional 



• Formation flight aims to improve 
induced drag at fixed lift... 


Cm — 


Cj 

tv ARe 


• Span efficiency factor, e, is a 
natural choice 

e = 3 

TV ARC Di 


Initial Mesh (100k cells) 
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Flow Analysis: Cart3D 


Requirements: automatic control surface deflections, geometry re- 
meshing, and high-fidelity CFD solutions 


• Inviscid analysis package, 3D 
compressible Euler equations 

• Unstructured Cartesian cells with cut- 
cells at wall 

• Adjoint-Based Mesh Adaptation 

• Highly automated for parameter 
sweeps, geometry manipulation, and 
re- meshing 
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Flow Analysis: Objective Functional 



• Formation flight aims to improve 
induced drag at fixed lift... 


Cm — 


Cj 

tv ARe 


• Span efficiency factor, e, is a 

natural choice „ 

Cl ~0.5 


e = 


cl 

tv AR C £>i 


C D ~ 0.008 


Initial Mesh (100k cells) 
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Flow Analysis: Objective Functional 


• Formation flight aims to improve 
induced drag at fixed lift... 


Cdi — 


Cj 

tv ARe 


• Span efficiency factor, e, is a 

natural choice „ 

Cl ~0.5 


e = 


cl 


TV ARC Di 


C D ~ 0.008 


Expensive to converge due to 
drag sensitivity 
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Results Outline 



• Mesh convergence study 

• Incompressible Effects 

- roll-trim 

- Vortex position sensitivity 

• Compressible Effects 

- roll-trim 
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Flow Analysis: Objective Functional 


• Formation flight aims to improve 
induced drag at fixed lift... 


ARe 


TV 


• Span efficiency factor, e, is a 

natural choice „ 

Cl ~ 0.5 


e = 


cl 


TV AR C Di 


Cd- 0.008 


Final mesh (20 mil. cells) 


Expensive to converge due to 



drag sensitivity 
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Formation flight metric 



drag fraction = j?;P in formation 

-^out formation 

• Provides a measure for overall drag savings of the formation 

• Focus on 2-aircraft echelon formations... 


drag fraction = 


-Plead H~-Pt rail 

2Plead 


^rrrrm,rr.iiiT~rri'r.i-rr-nni)rmTrr'jnr~n-rr-rrnirmTrr'jn!'T.nrr-rrn)^.iTT‘rr!'T.nrr-rrn)^.Trr‘rripmTrr'jn)rr.iTr‘rri‘r.--rT-m)Tr.irr‘n‘riiii, 

Ex.: a drag fraction of 0.9 represents a 10% drag savings 
: in formation flight 
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Results: Vortex in empty domain 


Flow conditions 

• Mach = 0.5 

• CL (lead A/C) = 0.55 

• Objective Function: Integrated 
pressure at 5 spans from inflow 

Goal 

• Verify mesh convergence of 
vortex in absence of geometry 


impose vortex BC 




5 

spans 

X 

trail aircraft 
CFD domain 
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Results: Vortex in empty domain 



inflow 
5 spans 



Inflow 
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Results: Vortex in empty domain 


Flow conditions 

• Mach = 0.5 

• CL (lead A/C) = 0.55 

• Objective Function: Integrated 
pressure at 5 spans from inflow 

Goal 

• Verify mesh convergence of 
vortex in absence of geometry 


impose vortex BC 


pressure 

sensor 


5 

spans 

trail aircraft 
CFD domain 
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Results: Vortex in empty domain 



inflow 
5 spans 
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Results: Vortex in empty domain 


inflow 
5 spans 
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Results: Vortex in empty domain 



Mach = 0.5 

• Objective functional: integrated pressure 5 spans downstream 
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Results: Vortex in empty domain 


inflow 
5 spans 
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• un-trimmed (baseline) configuration contains no ailerons 

• 1 -aileron deflected; increases lift on out-of-vortex wing 

• 2-aileron deflected; trims in conventional manner 


Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


36 


36 


Results: Simple Wing 


Geometry 

• Sub-sonic straight wing w/ 
NACA 001 2 sections 

Goals 

• QuantifyA/erify 
incompressible formation flight 
drag savings 

• Roll-trim effects 

• Determine trail A/C position 
sensitivity 



Geometry 

Cl 

Moo 

AR SrefUt 2 ) 

NACA 0012 Wing 

0.55 

0.5 

8 8 
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Results: Simple Wing 





Initial Mesh: 100k cells 


Final Mesh: 20 Mil. cells 
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Results: Simple Wing 


i \ 
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Results: Simple Wing 




Initial Mesh: 100k cells 



Final Mesh: 20 Mil. cells 
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Results: Simple Wing 



Final Mesh: 20 Mil. cells 
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Results: Vortex Position Sensitivity 




*View behind wing facing upstream 
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Results: Vortex Position Sensitivity 
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Results: Vortex Position Sensitivity 



• CFD point 


- 0.2 - 0.1 0 t 0.1 0:2 0.3 

y/b 


0.1 


0 


p 

N 


- 0.1 


*View behind wing facing upstream 
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Results: Vortex Position Sensitivity 
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Results: Vortex Position Sensitivity 


un-trimmed, optimal span efficiency 
2.2 @ wingtip 

Optimal vortex location migrates 
inboard w/ trim 


* 2-aileron trimmed outperforms 

1 -aileron 





HO 


-0.2 -0.1 0 ^ o:i 0.2 0.3 

y/b 


0.1 
0 1 
- 0.1 
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Results: Vortex Position Sensitivity 


Optimal vortex location migrates 
inboard w/ trim 

2- aileron outperforms 1 -aileron trim 

3- 5% decrease in drag fraction for 
2-aileron 



“ 0:3 

y/b 

Trimmed 1 aileron 

0.1 



"ID 5 0.1 0.2 03 

y/b 

Trimmed 2 ailerons 


:(? 


0.1 


• • 

0.2 -0.1 0 ~ 0.1 0.2 0 3 

y/b 


- 0.1 



Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


■ 

drag fraction: 0.70 0.75 0.79 0.84 0.88 


Results: Vortex Position Sensitivity 



• Largest induced rolling moments 
occur at wingtip and 20% overlap 

• Near 10% overlap, no induced rolling 
moments, yet close to optimal drag 
savings 
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Aileron Def.(deg.): -5.50 -2.75 0.00 2.75 5.50 48 
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Results: Vortex Position Sensitivity 
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Results: Vortex Position Sensitivity 




-20% -10% 0% 10% 30% 



y/b 
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Results: Vortex Position Sensitivity 
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Results: Vortex Position Sensitivity 




-20% -10% 0% 10% 30% 



y/b 
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Results: Vortex Position Sensitivity 





- 20 % - 10 % 0 % 10 % 30 % 
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Results: CRM 




Final Mesh: 40 Mil. cells 
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Results: Common Research Model 


Geometry 

• Modified CRM from 4th AIAA 
Drag Prediction Workshop* 

Goals 

• Quantify compressible 
formation flight drag savings 

• Determine roll-trim effects on 
more realistic wing/body 
geometry 



Geometry 

C L AR 

Srefift 2 ) 

CRM 

0.5 0.83 9 

4130 


*Vassberg, J. C.,et al. Development of a Common Research Model for Applied CFD Validation 
Studies, AIAA Applied Aerodynamics Conference, AIAA 2008-6919, August 2008. 
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Results: CRM 



" Out-of-Formation 


CRM at CL=0.50 
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Results: CRM 


" Out-of-Formation 


CRM at CL=0.50 



Mach Number 
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Results: CRM 



" Out-of-Formation 
• Formation (un-trim) 


CRM at CL=0.50 
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Results: CRM 


" Out-of-Formation 


CRM at CL=0.50 
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Results: CRM 



" Out-of-Formation 

• Formation (un-trim) 

• Formation (trim) 


CRM at CL=0.50 
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Results: CRM 


CRM at CL=0.50 
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CRM, Mach = 0.83 
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x/c 


x/c 
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CRM, Mach = 0.83 
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CRM, In-Formation, trimmed 
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CRM, In-Formation, trimmed 



Open Issues 



• Complete benefits map for summer conference 

• Extend analysis to > 2 aircraft formations 

• Effects of heterogenous aircrafts in formation 

• Deliver analytic drag model for formation flight including trim/ 
compressibility for use in NAS models 
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Summary 


• Performed relevant verification for vortex boundary condition 


• Developed methodology for evaluating formation flight benefits for 
2-aircraft echelon formations 


• Developed quantitative benefit maps for trail aircraft positioning 


• Evaluated erosion of formation flight benefits as a result of trim and 
compressibility 
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Wake Propagation Model 



Augmented Betz Method * 

• Far-field conservation method 

• Experimental core-size data 

• LES vortex decay rates 


Governing Equations: 


continuity: 

duy du z 

~dy ch~ 


vorticity: 
d( d( d( 
m +Uy d~ y + Uz d'z- p 


(a\ 

\dy 2 dz 2 ) 


Assumptions 

• Viscous effects neglected 
during rollup 


(where f = ^ 
V dy 


OUy 

dz 


Time Invariants: 


• All vorticity from lead is axis- 
symmetrically rolled-up into 2 
vortices 

*Ning, Andrew, Aerodynamic Performance of Extended Formation 
Flight, Journal of Aircraft, Vol. 48, No. 3, May-June 2011 
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r„ = 


r* = 


iy = 


j(„A 

J y(dA 
J z(dA 
J (y 2 + z 2 )(dA 
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Results: CRM at M=0.83 



Span Efficiency 


Results: Simple Wing 



cl 

7 r AR C jji 
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